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Abstract: Using the Hubbard model on a ID lattice and the "memory 
function" approach, we have calculated the electrical conductivity of a ID 
system of correlated electrons on a lattice. We have determined the temper- 
ature and frequency profiles of the conductivity. Effects of changes of the 
band filling on the conductivity were also briefly discussed. The results were 
qualitatively compared with experimental data on the Bechgaard salts, and 
the agreement is satisfactory. 



Introduction 

Studies of correlated electron systems are important for a number of 
reasons. In materials of reduced dimensionality, such as the quasi one- 
dimensional (QID) organic conductors, correlation effects play a highly im- 
portant role [1]. High Tc superconductors [2] and the fullerenes [3] are exam- 
ples of two and three dimensional materials in which, it is widely assumed, 
correlation effects are of crucial importance. 

The aim of this paper is to present results of a calculation of the contri- 
bution of electron-electron scattering to the electrical conductivity of QID 
organic metals. Theoretical studies of one dimensional correlated electron 
systems are interesting in attempts to explain and predict experimental data 
on such systems, as well as a possible source of insight in understanding the 
mechanism of high T^ superconductivity. The calculations will refer to the 
normal (that is, non-superconducting) state of QID organic metals, and they 
will be qualitatively compared with experimental results on the Bechgaard 
salts. The general chemical formula of these materials is (TMTSF)2X, where 
(TMTSF)2 stands for bis-tetramethyl-tetraselenafulvalene, and X is an anyon 
[4], [5]. 

A few years after the synthesis of the Bechgaard salts, it was shown that 
experimental results on their electrical conductivity can not be described 
by the standard theory of transport processes in normal metals [7] This 
result is the experimental motive for the present calculation. Following this, 
and a theoretical prescription by Anderson [8] , QID organic metals are 
usually studied within the Hubbard model. Although this model is one of 
the simplest in solid state physics, it is still a subject of intensive study (see 
[9]- [13] for some examples). General reviews of the field of organic conductors 
can be found in [14]- [16]. 

All the calculations in this paper were performed by the "memory func- 
tion" method. The main equations needed for the calculation of the conduc- 
tivity within this method are briefly presented in the next section, which also 
contains a derivation of the expressions for the real and immaginary parts 
of the conductivity. Details of the calculation are presented in the third 
part, while the fourth section is devoted to a qualitative comparison with 
the experimental data. 

The calculations were performed using S. Wolfram's MATHEMATICA 
software package, version 1. 2 on a PC 486 with 16 MByte RAM at 133 MHz. 



The method 

In studies of transport phenomena in various kinds of physical systems, 
one frequently enounters the problem of having to formulate and solve an 
appropriate transport equation. The memory function method gives the 
possibility of obtaining the conductivity without having to solve transport 
equations, but, instead, by representing the susceptibility in terms of a suit- 
ably chosen memory function. Details of the method are available in the 
literature (such as [17], [18] and references given there). 

The calculation of the conductivity is founded on the following two equa- 
tions: 



/•CO 

Xab =< A; B:^=-i exp{izt){[A{t),B{0)])dt (1) 

Jo 

and 

a{uj)=t{ujl/47rz)[l-x{z)/xo] (2) 



In eq, (2) Up = ArmeC^ /rrie is the square of the plasma frequency, rZe, e and 
rrie are the electronic number density, charge and mass , while xo = ^e/^^e 
is the zero frequency limit of the dynamical susceptibility x(^)- 

The integral in eq. (1) is a general definition of the linear response of 
an operator A to a perturbing operator B, and it is an analytic function for 
all non-real frequency z [17]. A(t) is the Heisenberg representation of the 
operator A. Inserting A=B=[j, H], where j denotes the current operator and 
H the Hamiltonian into eq. (1), one gets a definition of the current-current 
correlation function. The chemical potential is a function of the band-filling 
n, the inverse temperature (3 and the hopping integral t. The form of this 
function for the ID Hubbard model on a lattice with lattice spacing s has 
recently been determined as [20]: 



^ = [{[3tf{ns - 1) I t |]/[1.1029 + .imA{l3tf + mU{[3tf] (3) 



Taking the lattice constant s=l, and in the special case of half-filling 
(n=l), it follows from eq. (3) that /i = 0, which is in agreement with standard 
results [19]. As a first approximation, it will be assumed throughout this 
paper that the electrons form a normal Fermi liquid (however, compare [21] 
and [22]). 

The Hubbard Hamiltonian has the following form 



H = -t Y,{ct,aCi+i,a + c^i,aQ,a) + UY, nj^Uj^i (4) 

l,cr j 

which can be represented as H = Hq -\- Hi, the first term being due to elec- 
tronic hopping, and the second one to the inter-electronic on-site repulsion. 

The current operator is given by 



J = -^tJ2i4,aCl + l,a - ct^l,aCl,a) (5) 
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The functions a{z) and x{^) can be expressed in complex form as:(T(z) = 
(Tij(z) + i(7j{z) and x(z) = Xr{.z) + ix/ (-2) -Assuming that z= zi + iz2 , and 
Z2 = azi with a ;^ 0, it follows from eq. (2) that 



47r(l + ia)zi xo 

It follows after some obvious algebraic manipulations that the real and 
immaginary components of the electrical conductivity are given by the ex- 
pressions 



47r(l + a^)zi Xo Xo 



and 



_ ^P r. XR + aXl . f^. 



In the special case a = 0,these two formulas reduce to 



^PXl .Q^ 

ctr = (9) 

AttziXo 



and 



a, = ^(l-^) (10) 

Anzi Xo 

The calculation of the conductivity is in principle straightforward, but tech- 
nically complicated. The difficulties arise due to the fact that after inserting 
all the "ingredients" into eq. (1), one gets an almost intractable expression. 

The calculations 

In order to calculate the conductivity, one has at first to determine 
the current-current correlation function. It can be obtained by inserting 
A = B = [j,H] into eq. (1). 



The evaluation of this commutator is simphfied to some extent by the 
decomposition of the Hubbard Hamiltonian indicated after eq. (4). It can 
be shown that [j, Hq] = 0, and that 



^ = [i. H] = [j, Hi] = -itf/^(q+^Q+i,<, + q+ i_^Q,^)(5i+ij - Sij)nj^^^ (11) 

where all the symbols have their usul meaning. Transition to k-space can 
be performed by the relations of the following general form [6] : 



^l,a 



N^'/'J2^M^kils)cl,. (12) 
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The symbol N denotes the number of lattice sites, s is the lattice constant 
and L=Ns is the length of the specimen. Using the fact that in the case of the 
ID Hubbard model e{k) = — 2tcos(A;s),and introducing the temporal evolu- 
tion by relations of the form Cfc(t) = exp{—ie{k)t)ck , one gets the following 
final expression for x(z) 

x(^)=Ep,5,M(32^(V[(l+exp(/9(-/i-2tcos((?))))(l+exp(/5(-^-2tcos(A;))))] 
-1/[(1 + exp(/5(-/i - 2tcos(p))))(l + exp(/5(-/i - 2tcos(g))))](f/t)2 
{azi + i{zi + 2t(cos(g) + cos(p) — cos(5') — cos(A;)))) 
(cos(p + 5f)/2)(cos((g + A;)/2)[cosh(5(-p) - 1]/ 

{N'^iiaz-iY + (zi + 2t(cos(g) + cos(p) - cos(^) - cos{k)))Y)) (13) 



The summations were limited to the first Brioullin zone, and the lattice 
constant s has been set equal to 1, which means that L = N. Calculating 
the sums in eq. (13) is a non-trivial problem. They have to be evaluated 
for a 7^ 0,because this condition is built-in into the definition of the function 
x{z)- As the frequency is a real physical quantity, in the final expression for 
the dynamical susceptibility the limit a ^ will be imposed. 



Performing the summations in eq. (13) by MATHEMATICA, one gets 
the following approximate expression for the dynamical susceptibility: 

X = (32i(-l + cosh(l))[(l + exp{p{-fi + 2t cos(l - Tr))))-^- 
(1 + exp(/5(-/i - 2t)))-^]{Ut/Nycos\{l - 2n)/2)[{z2+ 
i{zi + 2t(-2 - 2 cos(l - 7r))))/(4+ 



{zi + 2t(-2 - 2 cos(l - 7r)))^+ < 2267 > 



(14) 



where << 2267 >>denotes the number of omitted terms. This equation 
is obviously untractable because of its length, and it has to be truncated after 
a certain number of terms. Limiting eq. (14) to its first 32 terms, multiplying 
out the products and powers, and expressing the result as a sum , it follows 
that the real part of the dynamical susceptibility is 

XR= [128UH^ cos^iil-2n)/2)]/[il+expiPi-fi-2t))YN\zl + {z^+2bty] 

+ < 527 > (15) 

and b = 2(-2 - 2 cos(l - tt)). 

Taking the limit Z2 ^ (because the frequency is a real quantity), and 
limiting eq. (15) to its first 20 terms, it follows finally that the real compo- 
nent of the dynamical susceptibility can be expressed as follows: 

Xr{^) = E^ KJ{lo + htf + E, L^uj/{u + htf (16) 

The explicite form of the functions i^jand Lj in eq. (16) can be read-off 
from eq. (15) developed to any given number of terms. 

The immaginary part of the dynamical susceptibility ( denoted by xi) can 
be calculated as the following integral [31] : 



XiM = -2{cuo/7c)PJ,^XR{^)dco/iuj' - ^i) (17) 

Limiting the summations in eq. (16) to terms with i,j < 4 using eq. 
(17) and taking the coefficients K and L from eq. (15) , one arrives at the 
following expression for xi'- 

Xi = {2bt/7r){Ut/N^f[iOo/{uJo + 2ht){ujl - (26t)2)] [4.53316(1+ 

exp(/3(-^ - 2t)))-2+24. 6448(1 + exp(/3(-/i + 2t cos(l - 7r))))-2] + 

{2/7r)[uj,/{ujl-{2htn{m/N^f 

[42.49916(1 + exp(/3(-/i - 2t)))-'^ + 78.2557(1 + 

exp(/5(-^+2tcos(l-7r))))-2] (18 

Inserting this result into eq. (9), one gets the following expression for the 
electrical conductivity of QID organic metals: 
(^rM = u;pX/ /47rcc;oXo = 

= il/2xo)iujl/iT)[ui-ibm-\Ut/NY{A2A9916{l+eMPi-f^-m))-'+ 
78.2557(1 + exp(/5(-/i + 2tcos(l + Tr))))-^ + (6t/(cuo + bt)) 

[4.53316(1 + exp(/5(-/i - 2t)))-^ + 24.6448x 
(l+exp(/?(-/i+2tcos(l+7r))))-2]} (19) 

Discussion 

Equation (19) is the final result of this paper for the real part of the 
electrical conductivity of QID organic metals. In order to test its validity, it 
has to be compared with real experimental data on the Bechgaard salts. 

This comparison will be performed with respect to three different exper- 
imentally controlable variables: 

- the temperature, 

- doping, and 

- the frequency. 



Model parameters in eq. (19) were chosen as follows: N = 150;U = 4t 
■,up = 3U; Xo = 1/3 and Uq^.GU. These values were at first chosen by analogy 
with high temperature superconductors. In the course of the calculations 
it emerged that similar values were used for the Bechgaard salts by other 
authors [25]- [27]. The lower limit for ujq was determined by imposing the 
condition cr^j > on eq. (19). 

Changes of the conductivity with the doping can be studied by varying 
the chemical potential , which, among other factors depends on the band 
filling. The case n = 1 (i. e. , /i = 0) corresponds to a half filled band 
and , experimentally , to a "clean" material. Deviations of the filling from 
the value n = 1 describe the doping of th specimen. Positive deviations 
correspond to doping by electron donors, while the opposite case (negative 
deviations ) describes the doping of a specimen by electron acceptors. 

As a first test. eq. (19) was applied to the case of a half filled band. 
Inserting n = 1 into this expression and developing in t as a small pa- 
rameter, one gets that the conductivity is approximately given by a^ = 
10^'^ {up/uYf^ {4.56 + .Spt), which is close to zero for physically acceptable 
values of the input parameters. This example can be extended to the level 
of a general conclusion - that the half filled ID Hubbard model is an insula- 
tor, in line with known results such as [22] or [29]. Relating to experiments, 
this implies that weakly conducting phases of QID organic metals can be 
described by a ID Hubbard model whose band filling weakly deviates from 
1/2. 

The behaviour of the conductivity as a function of the temperature for 
various values of the band filling n was studied by fixing all the parameters of 
the system, and then applying eqs. (3) and (19). The aim of the calculation 
was not to fit in detail the experimental data on any particular Bechgaard 
salt, but, instead, to reproduce the general trends observed in resistivity 
data on these salts. Examples of experimental curves p{T) (where p denotes 
the resistivity and T the temperature) are shown on the following figures. 
Data on {TMTSF)2ReOi and {TMTSF)2FS03 come from the author's 
work [23], while Fig. 3 is taken from [30]. Numbers near the curves on 
Fig. 1 indicate values of pressure, (kilobars), at which the data were taken. 
Theoretical curves p{T) calculated by expressions (3) and (19) for various 
values of the input parameters are shown on figures 4. - 10. Values of the 
parameters are indicated on the figure captions, and p = 1 at T= 116 K. 
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A comparison of figures 1 and 2 with figs. 4. -7. shows that data on 
{TMTSF)2Re04 and {TMTSF)2FS03 correspond to values of the band 
fining between .7and .9. Data on the Bechgaard salts containing the anions 
CIO4 and PFq are best described by the theoretical curve with n = 1.2. Fig- 
ures 8 and 9 are, by their form, similar to the form of the experimental curves 
p{T) for the salts of the {TMTTF)2X type (fig. 2 ). The general conclusion 
is that the calculations described in this paper reproduce semi-quantitatively 
the general behaviour of the experimental resistivity versus temperature data 
for the Bechgaard salts. A better quantitative agreement could be obtained 
by making an algorithm in which all of the system parameters could be var- 
ied, and the resistivity calculated at each step. A comparison of figures 7 
and 8 illustrates how a small change of the system parameters can induce a 
drastic change in its conductivity, which a known experimental fact. 

In order to study the behaviour of the conductivity of a Q ID organic 
metal as a function of the frequency, eq. (19) can be re-formulated as 

aK(cuo) = A[ui-{bt)^]-'[Q+btZ{u;o+bt)-^] (20) 

The symbols A, Q and Z denote the non-frequency dependent functions 
in eq. (19). Developing the last expression in ujq up to second order, one gets 

a^(cuo) = -A{Q+Z){bty^+AZ{bty^^u;o-A{2Z-Q){bt)-^ijl+ (21) 

Fitting eq. (21) to the measured frequency profiles of the conductivity, 
such as [24] - [27], one could determine the functions A, Q, Z and t. Taking 
the limit ujq ^ in eq. (21), it follows that 

(Tn{uo = 0) = -A{Q + Z){bt)-^ (22) 

which is the approximate static limit of the conductivity. 

Fixing all the parameters, and then varying the band filling, gives the 
possibility of investigating the effects of doping on the resistivity. Examples 
of changes of resistivity with doping are presented on figures 11. - 14. Values 
of input parameters are indicated on each of the figures. 
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Clearly, small variations of input parameters lead to large changes of the 
conductivity. Note that in all the examples shown on the figures the resis- 
tivity apparently tends towards a " saturation" value for band fillings greater 
than 1 (that is, in the case of doping a specimen with electron donors). 

On the other hand, for some values of the input parameters, doping with 
electron acceptors (i. e. , reducing the band filling) leads to a decrease of 
the resistivity. This imples that in certain regions of the j3 — t — p space, one 
could achieve increases of conductivity (or, perhaps, a transition to supercon- 
ductivity) by doping the specimen with suitably chosen impurities. Details 
of this problem will be discussed elsewhere. 

Conclusions 

In this paper we have performed a calculation of the electrical conductiv- 
ity of QID organic metals. The calculation was performed by the "memory 
function" method, using the ID Hubbard model on a lattice, and the results 
were qualitatively compared with experimental data on the Bechgaard salts. 
The correlation functions were calculated by definition, which is a distinct 
advantage over some previous work (such as [28]). Theoretical tempera- 
ture and frequency profiles of the conductivity are in good semi-quantitative 
agreement with experimental data. Note that the agreement is achieved when 
the band filling deviates by about 20% from 1/2, confirming earlier results 
[29]. The calculation was performed within the Fermi liquid picture, which 
is justifiable by the energy range we were interested in (compare [26], [30] ). 

Work discussed in this paper will be continued in the future along two 
different lines. On the purely calculational side, efforts will be made to take 
into account more terms in eq. (14). In improving the physics, a transition 
from the Fermi liquid to the Luttinger liquid description will be made. 
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